Tetraquark resonances with the triple flip-flop potential, 
decays in the cherry in a broken glass approximation 
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We develop a unitarized formalism to study tetraquarks using the triple flip-flop potential, which 
includes two meson-meson potentials and the tetraquark four-body potential. This can be related to 
the Jaffe-Wilczek and to the Karliner-Lipkin tetraquark models, where we also consider the possible 
open channels, since the four quarks and antiquarks may at any time escape to a pair of mesons. 
Here we study a simplified two- variable toy model and explore the analogy with a cherry in a glass, 
but a broken one where the cherry may escape from. It is quite interesting to have our system 
confined or compact in one variable and infinite in the other variable. In this framework we solve 
the two-variable Schrodinger equation in configuration space. With the finite difference method, 
we compute the spectrum, we search for localized states and we attempt to compute phase shifts. 
We then apply the outgoing spherical wave method to compute in detail the phase shifts and and 
to determine the decay widths. We explore the model in the equal mass case, and we find narrow 
resonances. In particular the existence of two commuting angular momenta is responsible for our 
small decay widths. 
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A long standing problem of QCD is the one of the ex- 
istence of localized exotic states and the corresponding 
decay to the hadron-hadron continuum. This is a diffi- 
cult problem both experimentally, since exotics usually 
decay to several hadrons, and theoretically, since exotics 
couple to wavefunctions of at least two quarks and two 
antiquarks (notice that a gluon couples to a quark and an 
antiquark). Nevertheless there is no QCD theorem pre- 
venting the existence of exotics, say two-gluon glueballs, 
hybrids, tetraquarks, pentaquarks, three-gluon glueballs, 
hexaquarks, etc, and the scientific community continues 
to search for clear exotic candidates . 

In what concerns multiquarks, they may exist due to 
different possible mechanisms. The perspective of mul- 
tiquarks as possible molecules of hadrons in attractive 
channels, where attraction is for instance due to quark- 
antiquark annihilation, has already led to the compu- 
tation of decay widths, which turned out to be wide 
[l|, Q • Here we explore another multiquark perspective, 
related to the Jaffe-Wilczek model [343 where diquarks 
are bound directly by four-body confining potentials, pro- 
duced by confining fiux tubes or strings. In this per- 
spective, the main effort of the scientific community has 
been to search for bound states below the threshold for 
hadronic coupled channels, avoiding the computation of 
decay widths [H . Aparently, the absence of a potential 
barrier above threshold may again produce a very large 
decay width to any open channel. However Marek and 
Lipkin suggested that multiquarks with angular excita- 
tions may gain a centrifugal barrier, leading to narrower 
decay widths P. 

Here we explore the possible existence of localized 
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FIG. 1: Triple fiip-flop Potential potential. To the list of po- 
tentials to minimize including usually only two different me- 
son pair potentials, we join another potential, the tetraquark 
potential. 



tetraquarks, or of resonances, co-existing with the con- 
tinuum of open meson-meson channels. To address this 
problem, we utilize the triple fiip-flop potential. The 
flip-flop potential for the meson-meson interaction was 
developed [8l-[l2j. to solve the problem of the Van der 
Waals forces produced by the two-body confining poten- 
tials [l3l - [l6l [isl [25| . Confining two-body potentials with 
the SU(3) colour Casimir invariant \i-Xj suggested by the 
One-Gluon-Exchange potential, lead to a Van der Waals 
potential proportional to V'{r)/r times a polarization 
tensor. This would lead to an extremely large Van der 
Waals force between mesons, which clearly does not exist. 
Thus two-body confinement dominance is ruled out for 
multiquark systems. Traditionally, the flip-flop potential 
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considers that the potential is the one that minimizes the 
energy of the possible two different meson-meson config- 
urations. This solves the problem of the Van der Waals 
force. Here we upgrade the flip-flop potential, consid- 
ering a third possible conflguration, the tetraquark one, 
where the four constituents are linked by a connected 
string [U [23. The three conflgurations differ in the 
strings linking the quarks and antiquarks, this is illus- 
trated in Fig. [1] When the diquarks qq and qq have 
small distances, the tetraquark conflguration minimizes 
the string energy. When the quark-antiquark pairs qq 
and qq have small distances, the meson-meson conflgu- 
ration minimizes the string energy. Whereas tetraquark 
binding [2l|, [22| has been investigated with the triple flip- 
flop potential, here we address the resonance decay width 
whith the triple flip-flop potential, in particular we study 
localized states existing above the threshold of meson- 
meson channels opened for decay. 

For the qiq2q3qi system, there is evidence in Lattice 
QCD at least for static quarks, that the hamiltonian. 



where the potential is decomposed in. 



t^l234 =C + V, 
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in a constant term C, a short range screened two-body 
part, say a one gluon exchange including. 



Vshort — ^ ] • A 



and in a long range conflning part, 

Vconf = +0-imm(ri,r2,r3,r4) 



(3) 



(4) 



where a is the string tension and L^^in is the string Icnght 
minimizing the energy of the colour singlet 91(729394 sys- 
tem. Lmin is the minimum of the three possible string 
lengths, the tetraquark string Icnght 

-^Ti234 = ^1 / + ^2 / + r3 // + r4 II + riii , 
but also a meson-meson string length, 

ij\/i3J\/24 = + r24 , 
and the other meson-meson string length. 



LMiiAhs = '''14 + '''23 , 



(5) 
(6) 
(7) 



where j is the distance between the points i and j and 
/ and // are the two Fermat-TorriccUi-Stcincr points of 
the tetraquark. 

However the hamiltonian of Eq. ([1]) is technically quite 
hard to solve. Moreover we would like to capture the 
essence of a potential where some dimensions are com- 
pact and others are open. Thus for simplicity we dedicate 
this work to the case where, 



meson-meson, V ^ 2 r 



Q, 



tetraquark, V = r + \ Z p 



FIG. 2: Simplified Potential. Assuming that the two internal 
mesonic coordinates are identical, we only have two three- 
dimensional variables. Nevertheless we remain with one con- 
fined or compact variable and one open variable. 



• all the quarks and antiquarks have the same mass, 
although we assume that they are different particles 
to avoid having to take account of exchange effects, 

• the short range interaction and the constant term 
are neglected in order to single out the effects of 
the conflnemcnt. 

• the quarks are non-relativistic for simplicity, 

• and we also simplify the number of variables, con- 
structing a toy model, in order to identify the 
physical mechanism possibly leading to localized 
tetraquark states. 

In Section II we detail our toy model, analogous to a 
quantum cherry in a broken glass, since our system is 
open to the continuum in some direction, and is conflncd 
in the other direction. In Section HI, we solve the cor- 
responding Schrodinger equation, searching for localized 
states although our system has no potential barrier and 
is open to the continuum. In Section IV we estimate the 
decay width of the tetraquark. In Section V we conclude. 



II. A TOY MODEL FOR THE POTENTIAL, 
SIMILAR TO THE ONE OF A CHERRY IN A 
BROKEN GLASS 

Classically, a sliding cherry in a glass is equivalent, for 
small oscillations, to a two-dimensional harmonic oscilla- 
tor potential. However if the glass is partially broken, as 
in Fig. [3l the cherry may escape from the glass, and this 
is similar to our flip-flop potential. Note that in quan- 
tum mechanics, the text book example of a metastable 
system, or resonance, is the one of a particle decaying 
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cherry in a broken glass - side view 



Cherry in a broken gloss - top view 



FIG. 3: Cherry in a broken glass. Our simplified two-variable 
toy model is analogous to the classical mechanics textbook 
problem of a cherry in a glass, but a broken one where the 
cherry may escape from. Here we solve this model in quantum 
mechanics, addressing the decay widths of a system compact 
in one variable and open in the other. 




FIG. 4: Plot of our simplified fiip-flop potential, as a function 
of the two radial variables r (compact) and p (open) . 



through a potential barrier with a tunnel effect. However 
the cherry in a broken glass is not well known in quantum 
mechanics, it is not clear a priori that it leads to reso- 
nances since energetically it is completely open for the 
escape of the cherry through the hole in the glass. Nev- 
ertheless, since some classical orbits of the cherry may 
remain trapped in the glass, it is interesting to study a 
quantum system in that class of potentials. 

We now simplify the triple flip-flop potential, in order 
to reduce it to a two-dimension problem, similar to the 
one of a cherry in a broken glass. It is convenient for 
a first study of the tetraquark resonances, to reduce the 
number of variables. While a convenient set coordinates 
for the 91(729394 system is, 
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ri2,34 
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ri - r2 

rs - r4 
miTi -I- m2r2 



mi + 1712 ma -|- 7714 

miri + TO2r2 + msrs + vtlaTa 
nil + ^2 + ^3 + TTiA 



(8) 
(9) 

(10) 
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for simplicity we assume in our simplified toy model 
that the diquark intradistance and the anti-diquark in- 
tradistance are similar, thus we assume that = ^34- 
Then we are left with only two distances, the diquark- 
diantiquark distance r and the intermeson distance p. 
Moreover, if we consider that on average p is perpendic- 
ular to r, and considering the effect of the Fermat angle of 
27r/3 in the tetraquark, our fiip-flop potential then sim- 
plifies to two different cases, the meson-meson case and 



the tetraquark case, 

VMM{r,p) = cr(2r) , 
Vrir, p) = cr(r + Vip) 



(12) 
(13) 



The corresponding flip-flop potential in this two- 
dimensional variable system, is simply. 



VpFir, p) = min(VA/A/, T^t) 



(14) 



and we depict it in Fig. HI 

In what concerns the kinetic energy, we assume that 
the momenta of the two Jacobi coordinates Pi2 ^'iid P34 
are identical. Thus we finally get the Schrodinger equa- 
tion, 



2m 



{^l + ^l)+VFF{r,p) 



i^{r,p)=E^{r,p) (15) 



and the unitary solutions of this equation constitute the 
main point of our paper. 



III. STUDYING THE SPECTRUM OF THE 
SIMPLIFIED HAMILTONIAN IN A CLOSED 
BOX WITH THE FINITE DIFFERENCE 
METHOD 

In the Schrodinger Equation (jlSp . separating the ra- 
dial from the angular coordinates, we get for the radial 
component. 
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r dr^ 

2 + -'^-^ V' 
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FIG. 5: The first localized state in the simplified 2- 
dimensional flip-flop potential with C = —1.0. 



and we discretize the equation with an anisotropic spac- 
ing in r and p, 

ri = ia , i = 0, • • • Nr — 1 , 

Pj^jb , j^O,---Np-l, (17) 

utilizing finite differences for the kinetic energies, 

(f 1 

u{ri,Pj) ^ {u^J-l - 2u,j +u,j+i) (18) 

where Uij = riPjipij. Thus the 2-dimcnsional 
Schrodinger equation is equivalent to a sparse matrix 
eigenvalue equation of dimension Nr Np x Nr Np. 

Notice that the boundary three-dimensional conditions 
for the radial equations are so that u{r, p) must vanish at 
the origin of both r and p, for the wavefunction ip to be 
regular. It is also convenient to utilize a box quantiza- 
tion, so the wavefunction also vanishes at the maximum 
distances an and bm wc reach, respectively for r and p. 

We solve the corresponding eigenvalue equation for a 
sparse matrix. In what concerns the spectrum, we get 
n m eigenstates of the hamiltonian. In the p direction, to 
have the best possible simulation of the continuum, we 
use a 6m as large as possible. 

Most of them are continuum- like states, i.e. they ex- 
tend to the end of the box, and in the limit of an infinite 
box their spectrum would be continuous, and only a few 
may be localized states. 

To search for the localized solutions we compute the 
Radius Mean Square (RMS) of each solution, in the open 
p direction. A localized state should have a smaller RMS 
than essentially all other states. Most states are expected 
to be continuum states with a RMS ~ 5m/2. In s- waves 
we find no localized solutions. But to check our formalism 
we can force the existence of a localized state by adding 
an ad hoc negative constant C to the tetraquark potential 
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FIG. 6: The first bound state in the simplified 2-dimensional 
flip-flop potential with C = —1.5. 




FIG. 7: Semi-localized state, or resonance for Ir = 1. 



aLTi23i ■ The corresponding partly localized state is de- 
picted in Fig. [5] While some continuum-like tail remains 
in the wave-function, indicating that this is a resonance 
state, a sizeable part of it is localized close to the ori- 
gin, where the tetraquark potential exists. If we further 
reduce the constant C we will get a confined state. In 
Fig. iniwe show a bound state we obtain for the constant 
C = —1.5 . Moreover it has been argued by Karliner and 
Lipkin Q that higher angular waves would create a cen- 
trifugal barrier that could originate resonances. Indeed 
we find resonances and even boundstates when we con- 
sider higher Ir- In Fig. [7]wc depict a resonance occurring 
a.t Ir = 1 and in Fig. [5]wc depict a localized boundstate 
occurring at = 3. 

To study the continuum-like solutions, it is convenient 
to measure their momenta, or wave number, k. With the 
momentum we can understand the spectrum and also 
measure phase shifts. To measure the momentum k and 
the phase shifts S in each of the non decaying channels, 
we have to measure the projection of the wavefunction 



FIG. 8: Bound state for Ir = 3. 



FIG. 9: Momenta of the various components as a function of 
the energy. 



in each one of the channels. Note that the wavefunction 
can be expanded as 



The arc the cigenfunctions of the hamiltonian 



Hmm = -— +2(77 
2m 



(19) 



(20) 



where we use equal quark masses, riii = m, and are given 
in terms of the Airy function 



Ai(- 



(21) 



where the are the zeroes of the Airy function. So, the 
ipi functions can be calculated by using the expression 



i/jiip) = / (fr(f>i{r)*^{r,p). 



(22) 



To measure the momenta kj and the phase shifts J.j, we 
simply fit the large p region of the non- vanishing ipi to 
the expression 



-0* -> A, 



sin(fcip + Si 



(23) 



As can be seen in Fig. [9l the momenta ki obey the 
relation 



(24) 



where is the threshold energy of the respective channel. 
Note that there are different thresholds opening at higher 
momenta, corresponding to radial excitations in the r 
variable. 

We also show the corresponding results for the phase 
shifts Si are shown in Fig. [TUl As can be seen, the be- 
haviour is irregular when we have more than one channel. 



this is due to the different contributions of multiple chan- 
nels, for each eigenstate calculated in the finite difference 
scheme. 

While the finite difference method is adequate to ob- 
serve the quasi-localized resonances and the localized 
boundstates, and also to measure the momenta of the 
continuum states, we now move on to a different method 
with the aim to study in detail the phase shifts. 



IV. STUDYING THE PHASE SHIFTS VIA THE 
OUTGOING SPHERICAL WAVE METHOD 

The irregularity of the phase shifts extracted from 
the finite difference solutions in a closed box, show in 
Fig. 1101 motivates us to utilize a precise and unitarized 
method to compute the phase shifts. We now solve the 
Schrodinger equation with spherical outgoing waves, to 
explicitly compute the phase shifts. 



A. method 

In a coupled channel problem, the Schrodinger equa- 
tion is. 



2m. 



(25) 



Considering the scattering from the channel i to the chan- 
nel j, wc have the following asymptotic relations, for the 
i channel, 



and, for j ^ i, 



(26) 



(27) 
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. n^ = 
. n'= 1 
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( E - 4 m c' ) / 
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in which case the equation becomes 



2m, dr^ 



+ V,ju = {E- e,)ui - V,hji{kr)r . (35) 



Thus, to calculate the phase shifts and the cross sections 
we only need to solve Eq. psp equation and to ana- 
lyze the assimptotic behaviour of each component of the 
wavefunction. 

We can reduce our problem in the dimensions p, r to a 
one-dimensional problem in p but with of coupled chan- 
nels. We just have to expand the two-dimensional wave- 
function as 



(36) 



FIG. 10: "Phase shifts" obtained from the finite difi^erences ( 
by projecting the eigenstates in the meson-meson eigenstates 
). As can be seen the behaviour is irregular when we have 
more than one channel, this is due to the different contribu- 
tions of multiple channels, for each eigenstate calculated in 
the finite difference scheme. 



if the channel j is open, otherwise vanishing. With the 
conservation of the probability, we obtain the optical the- 
orem, 



(28) 



This can be formulated in terms of partial waves, in which 
case the Eqs. and (P7|) become. 



u\ -> sin(fcir + In) + [21 + l)/-,e 



and 



4->(2Z + l)4e'' 
and the optical theorem becomes 



(29) 



(30) 



(31) 



The /j' can be computed by considering outgoing solu- 



(32) 



tions of the Eq. ([25]) , 

We have the equation for Xi 



- + V,jX3 ^{E- ti)x^ - V,he'^'^- , (33) 

and it can be simplified if Vij is spherically symmetric, 
writing Xi as, 



(34) 



where the 0^ are the eigcnfunctions of the r confined 
hamiltonian of Eq. (|20p . The one-dimensional potentials 
Vij are given by 

V,{p)^ j ^*{r){VFp{T,p)~VMM{m,{r) (37) 

where we subtract Hmm from the hamiltonian, since 
Hmm is already accounted for in its eigenvalues e^, to 
be used in Eq. (^5]) . 

Importantly, note that wc have two distinct angular 
momenta, which arc both conserved, = r x and 
Lp = p X Pp. So. each assimptotic state is indexed by 
its angular momentum Z,. and its radial number n^, and 
the scattering partial waves are indexed by I p. Thus the 
system can be diagonalized not only in the scattering an- 
gular momenta Lp but also on the confined angular mo- 
menta Lr. Wc can describe the scattering process with 
four quantum numbers: The scattering angular momen- 
tum Zp, the confined angular momentum 1^ and the initial 
and final states radial number in the confined coordinate 
r, rii and Uj. 



B. phase shifts 

We now compute the phase shifts, in order to search 
for resonances in our simplified fiip-flop model. Solving 
Eq. p5p for this system we can compute the partial cross 
sections and the total cross section for the partial wave I 
— either directly or by using the optical theorem of Eq. 
([51]) — and determine the phase shifts as well. 

Note that our flip-flop potential has the same scales 
of the simple Schrodinger equation for a linear potential, 
which becomes dimensionless when we substitute , 



/"ma\ 



1/3 



1/3 



E 



E 



(38) 



Thus the number of boundstates or resonances is inde- 
pendent both of the quark mass m and of the string con- 
stant a. Also, in the remaining figures, the energy will 
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FIG. 11: S-wave scattering cross sections from the channel 
with Ir = and rir = 0. 
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FIG. 12: S-wave scattering cross sections from the channel 
with Ir = 1 and rir = 0. 



be divided by the non-relativistic energy scale 



n =0 



(39) 



and thus our figures are dimensionless. Note that Eq 
is also the only energy scale we can construct with Ti, a 
and m. the three relevant constants in the non-relativistic 
region. 



1. the centrifugal barrier effect 

On Fig. [Tl]we can show the = partial cross sec- 
tions for the scattering from the channel with = 
and 71,, = 0. In Fig. [T^] we show the same results but 
for Ir = I. Interestingly, the bumps in the cross section 
seem to occur prior to the opening of a new channel. 

In Fig [13] we compare the phase shifts for different val- 
ues of Ir, namely for Ir ~ 0, 1,2 and 3. For Ir = 0, we 
don't observe a resonance, since the phase shift doesn't 
even cross 7r/2. However, for the Ir = 1 and Ir = 2 
cases, the phase shifts clearly cross the it/ 2 line, and a 
resonance is formed. This behaviour is somewhat ex- 
pected, since a centrifugal barrier in r would, in the case 
of a true tetraquark, maintain the two diquarks sepa- 
rated, favouring the formation of a bound state. The 
tendency of greater stability for greater orbital angular 
momenta seems to be further confirmed by the = 3, 
where besides the resonance, a true bound state seems to 
be formed, as can be seen by the different qualitative be- 
haviour of the phase shifts for this case. This bound state 
formation confirms our observation of a localized state, 
in Section III, with the finite difference simulation. The 
resulting wavefunction is presented in the Fig. |8l We 
also show the phase-shifts for scattering in the rir = 1 
and various Ir if Fig. 1141 ^'iid for scattering with Ir = 1 
and various channels in Fig. 1151 
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FIG. 13: Comparision of the phase shifts for Ir = 0, 1, 2 and 
3, with rir = 0. 
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FIG. 14: Comparision of the phase shifts for different Ir, with 
rir = 1. 
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FIG. 15: Comparision of the phase shifts for Ir = and 
different Ur- 



nl2 



— c = 





— c = 


-0.25 
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-0.75 
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-1.00 


— c = 


-1.50 




( E - 4 m c" ) / E„ 

FIG. 16: Effect of a different constant on the "tetraquark" 
sector of the potential in the phase shift. 



meson states, their non-rclativistic kinetic energy is, 



2. including an attractive constant 

As a verification of our method and of our numerical 
codes, we also study the effect of including an arbitrary 
constant in the tetraquark sector of the flip-flop potential. 
We consider a potential of the kind 



V = (7 min 



C). 



(40) 



The use of a negative value for C should result in favour- 
ing binding within the tetraquark sector and, possibly, 
the appearance of a resonance, or of a bound state, in 
the spectrum. 

In Fig. [iniwe see the results for the lp = phase shift 
for Ir — Q and = channel. As can be seen the phase 
shifts become greater for more negative values of C. For 
C = —0.25 we almost have the formation of a true reso- 
nance, with 5 passing 7r/2 but falling shortly after. How- 
ever, for C ~ —0.50 the formation of a resonance is clear, 
with a jump of tt in (5. If we further enhance the negative 
constant C in the potential, we then see the appearance 
of a bound state. The formation of a resonance and then 
of a bound state with the increase of the absolute value of 
C was already observed in the previous section, using the 
finite difference approach. The plots of the wavcfunctions 
of the resonance and of the bound state, are shown on 
Figs. [5]and|ni This expected result of increasing the at- 
traction, confirms that the similar states observed while 
increasing the angular momentum are indeed resonances 
or boundstates. 



3. a relativistic correction to the binding energy 

Another source for binding in present in relativistic 
corrections to the kinetic energy. A completely non- 
relativistic approach leads to an error in the kinetic en- 
ergy. Note that although we have different scattering 



E,{k) = 



fc2 
2^ 



(41) 



while the correct non-relativistic dispersion relation for 
the mesons should be 



E,{k) = Ahc' + — 



(42) 



where /i^ is the reduced mass of the system, depending on 
the mass of the mesons and thus depending on their bind- 
ing energy. In the full tetraquark problem this would be 
a little more complex, since we would have more degrees 
of freedom, namely Eij{k) = (Mi + Mj)c^ + with 

/iy = j^^'^j^j . . In our model we only need to consider 
the case Mj = Mi, since the two mesons are no longer 
independent. We introduce the relativistic correction in 
the mass 



M, 



where Mi is given by 



Mj = 2 m 



2c2 



(43) 



(44) 



and where Bi is the binding energy, given by the eigen- 
values of Hmm- The factor 1/2 comes from dividing it 
equally by the two mesons. So, to each state i corre- 
sponds the reduced mass 



4c2 



(45) 



Without this relativistic correction, our sistem would 
depend only in three constants Ti, m and a, and chang- 
ing the mass would only change the scale of the physical 
results, Eq. Now we have another constant, the speed 
of light, c, and this we will obtain qualitatively different 
results for different values of m/ \fa . 
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FIG. 17: Phase shifts, including the reduced mass correction 
for Ir = 0. 
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FIG. 19: Phase shifts, including the reduced mass correction 
for L = 2. 



1 = 1 



live phase shift. 
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(E-4mc')/E„ 

FIG. 18: Phase shifts, including the reduced mass correction 
for Ir — 1. 



Indeed we observe that the relativistic correction en- 
hances binding. In Figs. [T71 [TH] and [TH] for the Ip = 
phase shifts of the Ir = 1,2,3 and n,. = state as a 
function of the scattering energy, for different values of 
m/ y/a, are shown. When Ir = Owe see that for low quark 
masses, a resonance is present, but for a mass as low as 
2y/a it is already destroyed. For Ir = 1 we have a peren- 
nial resonance, which becomes thinner by increasing the 
quark mass. For Ir — 2 if we the mass is sufficiently low, 
we can observe, not only one resonance, but also a bound 
state and a second resonance. 



r 

2 



dS 
dE 



(46) 



The results without and with the reduced mass cor- 
rection are given in Tables H] and HIl Clearly, the width 
becomes larger for greater angular excitations, and for 
smaller quark masses. 

Notice that this is a toy model, but nevertheless we 
can estimate the physical scale of our decay widths. For 
instance in the case of light quarks, both our physical 
scales are similar m ~ ^/a ~ 0.4 GeV, Eq. [M] yields an 
energy scale of Eq — 0.4 GeV. Thus we get, for = 1, a 
decay width of the order 



Tl ~ 0.037 X 0.4[G'eV^] ~ 15MeV . 



(47) 



For tetraquarks composed of charm quarks, with a mass 
of approximately 1.5 GeV bottom quark, the energy scale 
is Eo ~ 0.26 GeV 



lOMeV , 



(48) 



and for a tetraquarks composed of bottom quarks, with 
a mass of approximately 5GeV, the energy scale is Eq ~ 
0.17 GeV and so the width is 



6MeV . 



(49) 



As for the case of = 2, the decay widths arc F ~ 
0.131 X 0.4 GeV ~ 52 McV for light quarks, F ~ 34 McV 
for the charm quark, and F ~ 23 MeV for the bottom 
quark. For = 3 the decay widths F are 140 Mev, 90 
Mev and 60 MeV respectively. 



4- computing the decay widths 



CONCLUSION AND OUTLOOK 



Finally we can compute the decay widths of the ob- 
served resonances utilizing the derivative of the respcc- 



We study tetraquarks with a simplified model, where 
the number of Jacobi variables is reduced. Nevertheless 
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(E 


- 4mc^)/Eo 


r / Eo 


1 




6.116 


0.037 


2 




6.855 


0.131 


3 




7.462 


0.352 



TABLE I: Decay widths as a function of Ir- 



Ir 




{E - 4mc'')/_Bo 


r / So 





1.0 


5.001 


0.039 




1.5 


5.096 


0.022 


1 


1.0 


5.659 


0.137 




2.0 


5.903 


0.093 




3.0 


5.990 


0.075 




4.0 


6.031 


0.064 




5.0 


6.053 


0.053 


2 


1.0 


6.194 


0.586 




2.0 


6.510 


0.259 




3.0 


6.634 


0.209 




5.0 


6.736 


0.171 




7.0 


6.777 


0.162 



TABLE II: Decay widths as a function of Ir, with reduced 
mass correction. 

our computations are fully unitary, and we are able to 
study resonances above threshold, including their wave- 
functions, the associated phase shifts, and their decay 
widths. 

We conclude that with higher orbital angular momen- 
tum it is easier to form resonances, thus the existence of 
tetraquarks is favoured when the orbital angular momen- 
tum is finite. This is consistent with Karliner and Lipkin 
Q who proposed an angular momentum mechanism for 
the binding of pentaquarks. Adapting their mechanism 
to tetraquarks, if we regard the tctraquark as a diquark- 
diantiquark system, the angular momentum generates a 
centrifugal barrier between the diquarks, impeding the 
recombination of the quarks with the antiquarks to form 
mesons, and thereby increasing the stability of the sys- 
tem. Moreover we find that the angular momentum of 
the confined coordinate r anticommutes with the angu- 



lar momentum of the continuum coordinate p, and this 
is crucial to prevent the fast decay to the continuum. 

When a relativistic correction is included in our model, 
we also find that the formation of resonances is favoured 
by low quark masses. This happens since lower quark 
masses produce a higher binding energy, with a higher 
meson to quark mass ratio favouring the localization the 
of resonances or boundstates. We indeed observe, for 
sufficiently large angular momentum {1^ ~ 3) and suffi- 
ciently small quark masses, that we have not only the 
appearance of resonances but also of bound states. 

Extrapolating the results obtained with this simplified 
potential to the similar triple flip-flop potential, and uti- 
lizing the typical scales of hadronic physics, we find plau- 
sible the existence of resonances in which the tetraquark 
component originated by a flip-fiop potential is the domi- 
nant one. The decay of a tetraquark to a meson pair turns 
out to lead to a narrow width, of the order of 10 MeV, 
when compared to the typical decay width of hadrons or 
the order of 100 MeV. 

For instance in the case of a S=2 and J=2 light 
tetraquark decaying to a pair of p mesons, the decay 
width computed here is sufficiently small to be neglected 
when compared with the effect of the rho meson de- 
cay width on the total decay width of the tetraquark. 
Thus this suggests that light tetraquarks exist, with de- 
cay widths comparable to the decay widths of normal 
mesonic resonances 

We also remark that the unitarized methods developed 
here are amenable to future studies of the full tetraquark 
problem with the triple flip-flop potential. 
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